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We propose a new mechanism which can lead to ferromagnetism in Hubbard models containing 
triangles with different on-site energies. It is based on an effective Hamiltonian that we derive in the 
strong coupling limit. Considering a one-dimensional realization of the model, we show that in the 
quarter-filled, insulating case the ground-state is actually ferromagnetic in a very large parameter 
range going from Tasaki's flat-band limit to the strong coupling limit of the effective Hamiltonian. 
This result has been obtained using a variety of analytical and numerical techniques. Finally, the 
same results are shown to apply away from quarter-filling, in the metallic case. 



I. INTRODUCTION 



Notwithstanding its long history the study of itiner- 
ant ferromagnetism is known as one of the most difficult 
theoretical problems. At the early stage the problem of 
ferromagnetism was studied within the mean field ap- 
proximation to the Coulomb interaction. However, in 
1960's Kanamori, § Gutzwiller f| and Hubbard @ in- 
dependently studied the effect of electron correlations on 
metallic ferromagnetism, which is left out in the mean 
field approximation. According to them correlation ef- 
fects are so important that itinerant ferromagnetism is 
possible only in restricted situations: a large density of 
states at the Fermi energy is not sufficient to realize fer- 
romagnetism. Nagaoka Q{ and Thouless || pointed out a 
different aspect of ferromagnetism. Studying the motion 
of a hole (and a few holes) in strongly correlated elec- 
trons, Nagaoka and other researchers noticed that the 
possibility of ferromagnetism is a delicate issue depend- 
ing on the lattice structure and the density of electrons. 



Quite recently several exact results have been found 
for ferromagnetism in Hubbard and related models; ||||] 
However those results are restricted to a special density 
of electrons. In this connection we note in particular 
that Tasaki has taken up the problem of ferromagnetism 
in a series of interesting papers fic|-|l3|| , and has studied 
several new models, in which the relevant band is flat 
or almost flat. For such models he succeeded in proving 
rigorously that the ground state is ferromagnetic when 
the electron density takes a particular value, namely 1/2 
(quarter-filling), in which case the system is insulating. 

The purpose of the present work is to reexamine the 
same problem from a slightly different viewpoint. Our 
approach can provide a new way to look at the problem 
of ferromagnetism and can thereby examine the stability 
of ferromagnetism, at least in some cases, without any 
restriction to the electron density. 

Restricting ourselves for a while to one-dimensional 



systems, we will consider a model which interpolates be- 
tween systems with large density of states and Nagaoka- 
like ferromagnetism, and we focus our attention on the 
following Hamiltonian: 

n = s E(i c ^ + H - c -) - f E i^-iA+i" + H - c -) 
+t ' E H-^ + d I+i>> + H - c -] + £ o E n i+h° 



+^E 



(1) 



where s is assumed to be positive, while Eq and t are 
more or less arbitrary and restrictions will be made when 
needed. Index j takes integer values, and sites with in- 
teger index we call 'c' sites, while those with half-odd 
integer index 'd' sites. We consider cases where the num- 
ber of electrons N e is equal or less than the number of 
'c' sites (n < \, where n = N e /L and L is the total 
number of sites). The ferromagnetic state has total spin 

S = S max = 

This two band model, which is depicted in Fig. 0, has 
more degrees of freedom than the regular Hubbard model 
and it should be able to exhibit more phases. Some spe- 
cial cases has already been studied in literature. For t = 
and £o = £fb = t' 2 /s — 2s the dispersion of one of the 
bands becomes flat, and the ferromagnetism in the so 
called flat band models has been discussed in the works 
of Mielke JTJ] and Tasaki. fo) Moreover, for e = 0, 
U = U' and t = — s the model is equivalent to what 
Miiller-Hartmann studied in the low density limit. E ] 
Finally, for e = 0, U = U' = +oo, t = -s and t' < t it 
was discussed by Long and coworkers. p6| 

Let us also point out that by applying a particle- 
hole transformation to this model we get a Hamiltonian 
where the hopping amplitudes change sign, Eq becomes 
— So + U — U' , the Coulomb interactions are unchanged, 
and the density transforms as n — > 2 — n. 

There are several reasons to believe that the ferromag- 
netic ground state should be stable in Eq. (fil) for a rela- 
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tively wide range of parameters. 

First, the basic building block is a triangle with two 
'c' sites and one 'd' sites, and the triangles are con- 
nected to form a chain (one dimensional model), or to 
form some kind of two (see Fig. |2j) or three dimensional 
lattice. Taking this triangle and putting two electrons on 
it, the triplet ground state is favored over the singlet for 
a relatively large region of the hopping t' , Sq and inter- 
action U. Especially, if Eo — t' 2 / s — s, then the high spin 
state is realized for any finite value of U. 

Second, as we already mentioned, for e = £fb the 
lower band is flat, and it has been proved that the 
model is ferromagnetic for any U > at quarter filling 
(n = 1/2). A nearly flat band also means large density 
of states near the Fermi level. 

Third, anticipating the results of Sec. [II, we mention 
that for U — +00 and large Eq the usual second order pro- 
cess giving rise to effective antiferromagnetic coupling be- 
tween spins on neighboring sites in the case of the regular 
Hubbard model is suppressed. Instead we have a third or- 
der process with an effective exchange cx —st n je\^ where 
the sign can be adjusted by the hopping s and it is ferro- 
magnetic for s > 0. This is also valid away from quarter 
filling. The mechanism described is clearly of kinetic ori- 
gin and resembles Nagaoka-effect. 

These arguments raise the following questions on the 
model in Eq. (g): 

1) Does the high spin ground state survives if we put 
many triangles together? In other words, how robust is 
the Tasaki ferromagnetism against a perturbation which 
destroys the flatness of the band? And to what extent is 
the strong coupling perturbational argument valid? 

2) Are there any cases in Eq. (H), where metallic fer- 
romagnetism is realized? 



II. SUMMARY OF THE RESULTS 

To motivate the reader's curiosity, let us start by show- 
ing in Fig. ^ the phase diagram we got in the insulating, 
quarter-filled case. We see that ferromagnetism with a 
correct excitation spectrum turns out to be a robust fea- 
ture of this model and is realized in a fairly large param- 
eter region. It is not surprising if we recall that in some 
loose sense our model interpolates between the Nagaoka 
like ferromagnetism of kinetic origin in the large-U limit 
and flat band ferromagnetism with large density of states 
at the Fermi level, and that they both favor the fully 
polarized, ferromagnetic ground state. This is very en- 
couraging since it shows that ferromagnetism can show 
up for cases that are more realistic than the flat band 
model or Nagaoka effect. We emphasize that near flat 
band models are multi-band systems, with the impor- 
tant consequence that the interaction becomes momen- 
tum dependent ( and is not constant as in the one-band 



Hubbard models) and this effective long range interac- 
tion is essential to stabilize the ferromagnetic state for 
small values of U. Furthermore the role of the triangles 
is not restricted to generate effective ferromagnetic ex- 
change, but also prevents ferrimagnetism |l^] (not being 
a bipartite lattice) and its importance has been already 
recognized by Takahashi pl| and Shastry et al. . Un- 
fortunately, for the itinerant case the situation is not so 
clear. Although the phase diagram is similar to that of 
the insulating case, the boundary effects due to low di- 
mensionality of the model are enhanced and make the 
analysis difficult. 

Finally, let us remark that although we studied the 
one-dimensional model in details, we believe that simi- 
lar conclusions would hold for higher dimensional models 
as well. 

To reach the above conclusions we analyzed the model 
in Eq. ([!]) with different techniques in the appropriate 
limits: 

a. Limit of U = +00 and Eo — ► 00 , any density. In 
this limit we use a canonical transformation to derive a 
low energy effective Hamiltonian and we show that fer- 
romagnetic exchange is generated in the third order. 

b. Large U and Eo, quarter filling (n = 1/2): 
Rayleigh-Schrodinger perturbation theory allows us to 
derive a series expansion for the dispersion of the spin 
flipped state up to 7th order. The series is then analyzed 
with the help of Pade approximants and we determined 
the boundary of the local stability within the radius of 
convergence. In the lowest order the model is ferromag- 
netic for U >U C ^ sel/t' 2 . 

c. Small U and nearly flat band, quarter filling. For 
the cases when both U and the width of the lower band is 
small comparable to the band gap, we concentrate on the 
Hamiltonian restricted to states in the lower band. We 
find that there is a very nice Goldstone mode (magnon) 
with dispersion u>(q) = Dq 2 +0(q ) for q — > 0, where the 
spin stiffness D can be calculated analytically. Further- 
more, the ferromagnetic state is stable for U > U c oc \p\, 
where p measures the flatness of the band and is defined 
as p = e - £ F b- 

d. Quarter filling, arbitrary U and Eq, small size sys- 
tems. We diagonalize small clusters in the S z — S max — 1 
(spin flipped) subspace up to 20 sites and in the S z = 
subspace up to 16 sites. The data are consistent with the 
results obtained in the aforementioned limits, and we fill 
the phase diagram for any Eq. The global stability seems 
to indicate that the transition is from ferromagnetic to 
singlet state, except near the flat band where we do not 
have reliable data due to convergence problems. We find 
also that near eo = — 2s the system is very sensitive to 
boundary conditions due to orbital degeneracy of the fer- 
romagnetic state. 

e. U = +00 and large Eo, less than quarter filling 
(n < 1/2) Using the factorized wave function, we de- 
termine the sign of the effective spin-spin interaction for 
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different values of s, t' and t. We discuss the excitation 
spectrum of the spin flipped state for J c g <C s, t', t, and 
show that in case of orbital degeneracy of the ferromag- 
netic state we need a finite value of J e ff to make the high 
spin state the lowest in energy. 

/. Arbitrary U and Eq, eighth filling (n = 1/4,), 
small size systems. From exact diagonalization of L = 
12, 16, 20, 24 site systems in the S z = S max — 1 subspace 
with N e = L/4 electrons we determine the boundary of 
the local stability. Due to large finite size effect the re- 
sults are not always conclusive. The effect of orbital de- 
generacy is enhanced compared to the insulating quarter 
filled system. 



III. LOW-ENERGY EFFECTIVE HAMILTONIAN 
IN THE STRONG COUPLING LIMIT 

In this section, we will derive an effective Hamilto- 
nian when t, t', s <C £o, U, U'. In this limit, we have to 
distinguish between two cases depending on the relative 
magnitude of e to U and U': £q < U, U' and U, U' < £o- 
First we study the former case and we take the large- 
U limit, i.e. U,U' — > +oo. We mainly discuss the case 
where the electron density satisfies n < \. Since Eq is 
positive, only the lower band consisting mostly of the 'c' 
component is filled. 



For U, U' — * oo the Hamiltonian is reduced to 

n = s E + H - c -) - t E + H 



3,° 



(2) 



where c 9CT and d i are annihilation operators under the 

constraint of no double occupancy: c,j a = Cj a (l — rija), 
d i = d i (1 — n i _) and a = —a. The creation 

3-2° 3-2 a 3~2 S 

operators are defined similarly. 

To derive an effective Hamiltonian to describe low- 
energy physics, we can treat the t' <C £o term as a pertur- 
bation inducing transitions between the lower and upper 
bands. We use the canonical transformation to eliminate 
the off-diagonal processes due to 

H' = t'Y\(d ] i +d f i )£,-„+ H.c.1. (3) 

The other terms in Eq. (||) are denoted as TCq. Since 
the canonical transformation is standard, we skip the de- 
tails. Eliminating the first-order effect of t', an effective 
Hamiltonian up to the second order in t' is obtained as 
follows: 



(m'\H e &\m) = (m\H \m)6 mm ' + - 



(to' \H'\n){n\H'\ to) (to' \H'\n) (n \H' | to) 



E„ 



E n 



Em' E n 



(4) 



where to, to', n are eigenstates of H.q. The energy denominator E m — E n can be written as E m — E n = —Eq + E' rn — E' n . 
Here E' m represents the part due to t and s in the eigenvalue E m , i.e. H' \m) = E' m \m) where H' denotes Ho except 
for the £o term. Since £o is assumed to be much larger than t and s, we can expand the energy denominators in the 
second term of Eq. (^) in terms of H' Q as 



(m'\H eS \m) = (m\H Q \m)S mm , ~ —{w!\H!V!\m) + {m'\W [H ,H'} \ 

E IEq \ 

-X((m'\H'[H , [H' ,H'}}\m) + (m'\[[H' ,H ],H ]H'\m) 

Z£ \ 



{m'\[H',H' ]H'\m) 



(•5) 



As mentioned before, we assume n < |. Then only 
the lower band is filled. Having such states in mind, we 
evaluate each term in Eq. (ph. The result is as follows: 



3" 



zndH$=H% ) +Hi 3 s b) with 



o/(3a) _ _ 
^ •cff■ — 



2/7 



(6) 
(7) 



(3b) 
iff — 



St" 



+(%'-1<t + n j+la )(l - n js ) 



(8) 



(9) 



Although it looks complicated, Eqs. become 
simple for n = \, in which the states |to) are written as 
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TO) 



IR 



(10) 



where the product is taken over all integer sites and where 
the quantum numbers m just correspond to the set of 
spin indices cr(J). In fact it is easy to see that except for 
constant terms 7i c ff is reduced to 



K 



cff 



-4s 



3 



(11) 



where S,- 



jj means a spin | at site j. Eq. (^Tj) shows 
that the effective interaction is third order in the hop- 
ping; it is essentially a ferromagnetic Heisenberg model 
for positive s. Therefore the ground state is ferromag- 
netic with the maximum total spin. The reason for this 
result is simple. The contribution in Eq. (^Tj) comes from 
the ring exchange among 3 sites (a triangle formed with 
j, j + l,j + I sites, see Fig. |j) with the hopping having 
the same sign, which is known to lead to a ferromagnetic 
exchange. Needless to say the system is insulating with a 
charge gap of ~ e . Here we note that although Eq. ( p"l] ) 
has been derived for the one-dimensional case, it is appli- 
cable to more general cases like the ones shown in Fig. [|. 
In that case the summation over j should be interpreted 
as the one for all the nearest neighbor 'c' sites. 

So far we have assumed U, U' — > oo. Relaxing this con- 
dition, let us take into account the lowest term in 1/U. 
This is essentially the same as the one-band Hubbard 
model in the large-?/ limit. [^0) The effective Hamilto- 
nian which should be added to Eqs. is 

4s 2 x 

^off = — (Si • S i+ i - \nin i+ i) 



S 

"u 

1.(7 



H.c. 



(12) 



where n. 



is the density of electrons at site j, 



and for quarter filling it simplifies to 



^eff = — 



U ^ 



(Si • Si+i 4- 



(13) 



This Hamiltonian favors antiferromagnetic fluctuations 
and works against ferromagnetism, so the final state will 
be decided by the competition between this Hamiltanion 
and that of Eq. (§). 



IV. STABILITY OF FERROMAGNETISM AT 
QUARTER-FILLING (n = 1/2) 

We now examine how robust ferromagnetism is at 
quarter-filling. We will look mainly at the stability 



against one spin flip; this is the local stability condition. 
At some points we will mention the global stability as 
well. Henceforth we assume, for simplicity, U' — U and 
t = in Eq. (|l|) , in which case the topology of the model 
is the A chain. 



A. Large U and eo 

As we noted in the previous section, for eo positive, 
and for U and eo much larger than the hopping ampli- 
tudes, the quarter filled system exhibits a charge gap pl| ] 
A c min(C7, eo). Every 'c' site is occupied by one elec- 
tron, the excitations in the charge sector are pushed up 
to high energies, and the only low lying excitations are 
coming from the spin degrees of freedom, which can be 
described in lowest order by the competition of Hamilto- 
nians (|ll|) and |[|). 

To get reliable results in this limit, it is desirable to 
extend the calculation of the effective spin Hamiltonian 
to higher order corrections. Although in principle it can 
be done, see for example in Ref. [Q, it is a very diffi- 
cult task, especially for a relatively complicated model as 
in Eq. (Q). Therefore here we follow a slightly different 
approach, and we compute the dispersion relation of the 
spin flipped state given by 



s(k) = -Ki [1 - cos(fc)] - K 2 [1 - cos(2fc)] 



(14) 



where the coefficients are related to the exchange inte- 
grals in the spin Hamiltonians. Were the spin Hamilto- 
nian of the form J\ J2i Si ■ Si + i + J2 J2i Si • Si+2 + ■ • ■ then 
A'i = Ji, K2 — J2 would hold. However, the canonical 
transformation generates terms [E2| like (SjSj/)(SjSj/) 
etc., and this makes the correspondence difficult. 

The series expansion for the coefficients in Eq. (|l4| ) 
can be obtained using the standard Rayleigh-Schrodinger 
perturbation theory. As an initial state one picks up a 
state where all the 'c' sites are occupied with spin up 
electrons except for one site, where the spin points down- 
wards, furthermore all the 'd' sites are empty. Since the 
system is translationally invariant, momentum is a good 
quantum number, and the problem is simplified as the 
Hilbert space is spanned by states with given momen- 
tum k. In this way we avoid the L/2 fold degeneracy of 
the ground state, as in every k subspace the lowest en- 
ergy state is unique. Now, through a number of virtual 
processes (hopping to 'd' sites, or hopping to already oc- 
cupied sites), the spin down electron will finally return 
to a 'c' site. The distance j which the down spin elec- 
tron has traveled will pick up a phase factor e lk: > and will 
contribute to Kj. Although it does not look difficult, 
the technical realization of this calculation for higher or- 
ders is not straightforward. We used a variation of the 
method given by Barber and Duxbury, [ p3[ which is rela- 
tively simple. Recently, a more efficient method based on 
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linked cluster expansion become available due to Gelfand. 
11 

We have calculated the coefficients K\,K2, ... up to 
7th order, and they are given in the Appendix. The low- 
est order contributions are 

K w If! 



(3) 



4st'- 



8st' 2 
Ue ' 



K 



in agreement with Eqs. ( |ll| ) and (|13|). 

The first contribution to Ki comes in 4th order, 



(15) 



K, 



(4) 



= - 



As*_ 
773 ' 
16sH' 2 8s 3 t' 2 



e 2 U 2 



(16) 



and K3 appears first in 6th order (see Appendix for de- 
tails). 

To be more specific, let us now turn to the special case 
of t' = s. The series can be analyzed with Pade approxi- 
mants. A typical plot of K\ as a function of s/eq is given 
in Fig. ||. Fixing the ratio U /sq, we look for which value 
of Eq the coefficient K\ changes sign. It is enough to 
consider K\ since \Ki \ <C \K\\ holds near the transition. 
Changing the ratio of U/sq, we are able to determine 
the region where ferromagnetism is locally stable. The 
result is shown in Fig. ^. For s ^= t' there are some quan- 
titative changes in the location of the boundary, which 

(2) (3) 

can be already recognized from K \ and K\ yielding 
U c ^se 2 /t' 2 . 

A similar calculation can be carried over for U, —Eq 3> 
s, t' . Then the electrons sit on 'd' sites and for t = the 
first contributions to K\ come in fourth order: 



21 1 ~ -2tt 



K 



(5) 



4^ 
8t 5 



Hi 1 



m 5 



uy 



m 5 



]U el{2e + UY 



(17) 



where we concentrated on s = t' = t case, £0 = — £0 an d 
the coefficients are denoted Kj. For higher order correc- 
tion up to 7th order we refer the reader to the Appendix. 
The series is convergent for |eo|/i <^ 15. 



B. In the vicinity of the flat band 

Our model is equivalent to Tasaki's flat band model 
p0| for e = £fb = t' 2 /s — 2s, where he has rigorously 
proven that the model is ferromagnetic for arbitrarily 
small value of interaction. In this section we restrict 
£0 to be close to £fb, and we introduce a quantity p 
which measures the perturbation from the flat band as 



p — £0 — £fb- In a series of recent papers M Tasaki 
extended the existence of ferromagnetism for such per- 
turbed models, reaching the result that ferromagnetism is 
stable for the nearly flat bands as well for U > U c , where 
U c is roughly proportional to the perturbation strength, 
in our case p. His treatment is rather involved and we 
are looking for a simpler description. 

Let us consider the model Hamiltonian in k space: 

H = 2S C S kc L C ka + £ oYl d i° d k* 



2t'J2^s-(c{ a d ka +dl r7 c ka 



2U 



X! ( C L T c L 1 c k 3 1 c k 4 T + <4i T d k 2 1 d k 3 1 d k 4 T 



/Cl,/C2,&3,fc4 



(18) 



where in the summation fci + fc2~ — &4 = G should hold, 
G = 0, ±2ir is a reciprocical lattice vector. In the non- 
interacting case (U = 0) it describes two bands where 
the lowest band is nearly flat (the bandwidth is cx p), 
while the upper band is dispersive. The two bands are 
separated by a band gap Aband = 2s + £ + O(p)- We 
expect, as already noted by Kusakabe and Aoki, |13| that 
for small values of U the essential physics is now going 
on in the lower band. To construct an effective Hamilto- 
nian, we first diagonalize the hopping part of the Hamil- 
tonian by the following canonical transformation (here 
we choose t' > for convenience, as the model does not 
depend on the sign of t'): 



C L= cosa fc 4 CT +sino! fe &|. CT , 
d L = -sinatfcok + cosa fc 6^, 



(19) 



where a\ a and b\ a are the creation operators of elec- 
trons on the lower and upper band, respectively, is 
determined via 



tan 2ai 



At' cos(fc/2) 
£0 — 2s cos(fc) 



(20) 



and is a continuous function of the momentum k, so that 
it is between and ir/2 for momenta in the Brillouin 
zone; furthermore a±„ — for £0 > —2s and a± n = § 
for £0 < —2s. It is convenient to extend its definition 
outside the Brillouin zone so that «27r+fc = —oik for any 
k in case of £0 > —2s and a-i^+k = it — for £0 < —2s 
. The band dispersion is then given by 



£ a ,b{k) 



s cos k =F 



£0 
2 



s cos k 



At' 2 ( 



' 2' 
(21) 



the effective Hamiltonian in the lower band reads 
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H = ^e a {k)al cy a kcJ 

k.a 

+ E v ( k i, k 2,h,k 4 )al^al 2l a k3l a k4h (22) 

where the interaction is defined as 
2U ( 

V(ki,k 2 , k 3 , k 4 ) = —(cosa kl cosa k2 cosa ks cosa ki 

+e tG / 2 sina^j sinafe 2 sina^ sina^ (23) 

and ki + k 2 — fc 3 — k 4 = G holds. The factor e lG l 2 is 
1 for normal scattering (G = 0) and —1 for umklapp 
scattering (G = 2ir), and it can be conveniently taken 
care of using the definition of a k for k values outside the 
Brillouin zone. 

Here we neglected terms which are coming from the 
scattering on the states in the upper band [|TJ and they 
are of the order of J7 2 /Ab an d, so that the approximation 
is valid for U -C Aband ■ 

We now concentrate on the excitation spectrum of the 
spin flipped states of momentum q which are defined as 



lx(9^))=E/p^)4 +9 ,l a P ,Tl FP > 



(24) 



where |FP) is the fully polarized state with all the spins 
up (i.e. S — 5inax = L/2), q is measured from the mo- 
mentum of |FP) and v is some quantum number. Since 
the lower band is filled, these states span the whole re- 
duced Hilbert space and no approximation is made in 
solving the effective Hamiltonian (|22]). 

The dispersion of the spin flipped state uj(q) is given 
by the eigenvalue problem |13[ 



(25) 



where the matrix is given by 

Mp, p <(q) = [e a (p + q) + A(p + q)~ e a (p)} 6 P , P > 

-V(p,p' + q,p + q,p'). (26) 

A(p) = J2 k V{k,P,P, k) is the Stoner gap. The size of the 
matrix is {L/2) x {L/2) and can be easily diagonalized 
numerically for large system sizes. 

A typical plot of the excitation spectra is given in 
Figs. ^ and ||. There we can see that a low energy 
magnon (Goldstone mode) and an optical magnon (with 
spin oscillation within the cell) emerge from the Stoner 
continuum. For small values of q the magnon has a dis- 
persion io {q) = Dq 2 . 

Now, what are the conditions for local stability? First 
of all, to have a dispersion relation appropriate for ferro- 
magnets, stiffness D should be positive. The interaction 
strength where D changes sign we call Ud- However, this 
is not enough, as it can happen that although D > 0, the 



Stoner continuum pushes down the spin wave mode at 
q = it. Therefore the decisive criterion for local stability 
is that all the excitation energies are positive (i.e. all 
the energies are above the energy of the fully polarized 
state), and this will define U c . 



1. Calculation of the Ud 

Clearly D is given as \d 2 uo {q) / dq 2 \ q= o and it can be 
calculated from Eq. ( p6| ) by using Hellman-Feynman rule 
(it will give D-y) and second order perturbation theory 
{D 2 ), so that D = Di + D 2 , where 



D 1 = l{x(oM^P-\xM) 



9=0 



Do 



E^Kx(o^)l^lx(o,o))| 2 



lo v — UJq 1 



9=0 



(27) 



Here |x(0, 0)) is the ground state wave function for mo- 
mentum q = 0, which is nothing else but S^IFP) and is 
trivially given by /°(0) = \J2/L. The |x(0, v)) are ex- 
cited states. Since by definition D 2 is negative and works 
against ferromagnetism, sufficiently large D\ is needed to 
compensate D 2 and to give D > 0. 

Nonzero contribution to Ui will come from the off- 
diagonal matrix elements of M PtP i : 



i r dp r dp' d 2 

2j_^]_T T 2^dq~ 



2 V{p, P ' 



q,p + q,p 



9=0 

(28) 



since the diagonal part vanishes being an integral of the 
derivative of periodic function. Using the special form 
of interaction like in Eq. (^) and carrying out partial 
integration, one can show that this quantity is positive, 
unless V is a constant. In other words the momentum 
dependence of the interaction is crucial to stabilize the 
ferromagnetic state. 

Let us now turn to D 2 . In our case the nonvan- 
ishing matrix elements come from odd-momentum ex- 
cited states labeled by momentum p, where /|(0) = 

(<5p,p — 3 P ,- P )/"v2 with energies w v — cjo = Thus 
we obtain 



D 2 = - 



12y, 1 

4L^ A(p) 
p 

i r dp i 
~l]_^Mp) 



d 



de{p) , dA{py 2 



9=0 



dp 



dp 



(29) 



The divergence of D 2 for small values of U is the signa- 
ture that the dispersion becomes linear in q for [7 = 
and the result is non-perturbative in U. 

To make the calculation simple we expand everything 
in powers of the small parameter p: 
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t' 



COS Oifc 



smoj 



Vi' 2 + 4s 2 cos 2 (fc/2) 
2scos(fc/2) 



O(p), 
O(p), 



yjt 12 + 4s 2 cos 2 (fc/2) 
e(fc) = -2 + (1 - cos 2 a fe )p + C(p 2 ). 

The following integrals are useful for our purpose: 

r dk 2 f 

2tt Vi' 2 + 4s 2 

dk ■ 2 , r 

4 = / —Sin Qjt = l- J c , 



(30) 



^ dk f dcosctk 



2tt V dfc 



2i'(i' 2 +4s 2 ) 3 / 2 ' 



2?r V rffc 



2tVi' 2 + 4s 2 



(31) 



With the help of these integrals, we can write the Stoner 
gap as 

A(fc) = U(l c cos 2 a k + I s sin 2 a k ) . (32) 
Finally, stiffness D is the sum of the D 1 and D 2 - 



D 1 = U(I C I' C + IJ' S ), 

D 2 = -^j[-2p+{2I c -\)U] 2 , 



where 



C = 



dk 



cos afc 



ff 27r I c cos 2 a*; + I s sin a k 



d COS Qffc 

dk 



(33) 



(34) 



Although the coefficient C can be given in a closed form, 
it is rather complicated and we give its value for some se- 
lected values of t'/s in Table. |, together with the values 
of U where D changes sign (Ud)- Note that the equation 
we solve is a quadratic one and it gives two solutions: 
one is for p > and the other one is for p < 0. 
In the limit t' > s we get 



(35) 



D -K*-4 (f) 



for small U and p, which is in agreement with Tasaki's 
result. Ol 



2. Calculation of U c 

From numerical diagonalization of the matrix M p y we 
see that the other local minimum of the acoustic magnon 
branch is at q = n. Unfortunately we do not know how 
to get U c , therefore we have to solve the eigenvalue prob- 
lem of Eq. (|2^) numerically and check for which values of 
U u}(q = 7r) > holds. From this we obtain U c . Again, 



we have to distinguish between p > and p < 0. The 
results are presented in Table. |. 

We also learn that for small t' a very good estimate 
of U c can be obtained by using only the diagonal ele- 
ments of matrix M PjP /(7r) applying the Stoner criterion. 
Namely, the spin down electron occupies states in an ef- 
fective band £j.(p) = £ a {p) + due to the background 
of spin up electrons. This band is higher in energy than 
the £|(p) = £ a (p) of the spin up electrons. Accord- 
ing to Stoner's theory, when lowering U the Stoner gap 
decreases, and at some point we start to fill the band 
of spin down electrons. If we neglect the off diagonal 
matrix elements describing the interaction between the 
hole in the spin up band and the spin down electron, 
this happens exactly when the energies of the two bands 
start to overlap. For p > the minimum of £j.(p) is 
at p — 7r and the maximum of £|(p) is at p = 0, so 
that the criterion is £i(tt) — £^(0). This will give us 
jjapp _ 4p/(iV<' 2 + 4s 2 ); for p < the corresponding 
equation is £j.(0) = e^ 71 ")- 

For t' ^> s, on the other hand, Ud and U c are very 
close to each other. This is due to the fact that the ex- 
cited states in this case are far from the acoustic magnon 
branch and their interference is small for values of U close 
to U c . This was also predicted by Tasaki jll| in the sense 
that lower and upper bound on the dispersion relation are 
very close to each other in this case. 

We have compared these results with the exact diago- 
nalization results for the stability of the spin flipped state 
in Fig. H for t' = s and we found that U c increases linearly 
with p even for relatively large values of p (w 0.5s)for 
p > 0, whereas for p < the linearity holds only for 
small values of p. 



C. The orbital degenaracy at eq 



-2s 



A very interesting point of our phase diagram is at 
£q = — 2s, where the upper and lower bands touch. Then, 
depending on the system size (L — 8, 12, 16, 20 ... with 
periodic boundary conditions or L = 6, 10, 14, 18 ... 
with antiperiodic boundary conditions), the fully polar- 
ized ground state can be degenerate. This orbital degen- 
eracy seems not to favor the highest spin states, as it 
was noticed e.g. in the case of Nagaoka ferromagnetism 
by Siito. |23] To show this explicitly, we calculated the 
energy dispersion in the S z = SVax — 1 subspace (see 
Fig. 10) and find that the difference between the energy 
E of the groundstate in the S'max — 1 sector and the en- 
ergy Efm of the ferromagnetic state satisfies the following 
finite size scaling: 



E-E FM = ^+0(l/L 4 ) 



(36) 



with C < 0, thus for any finite system the ground state 
has a total spin less than S'max- 



7 



It seems that this strange behavior is constrained to 
to = —2s and away from it the S max state is the ground 
state in the thermodynamic limit for U > U c , To clarify 
this issue more precisely further work is needed. 



V. FERROMAGNETISM IN THE METALLIC 
CASE (n < 1/2) 

Here we are mainly concerned with the case of large U 
and £o . A nice feature of the effective Hamiltonian is that 
as far as the one-dimensional case is concerned, one can 
also discuss the metallic case corresponding to n < 1/2. 
(A different way to show the stability of ferromagnetic 
state is given in Ref. p6| .) Actually the metallic case is 
more interesting from a physical point of view than the 
insulating state. For n < 1/2 all the terms in Eqs. (D)~(||) 
contribute. Although it looks complicated, the effective 
Hamiltonian in Eqs. (|])-(|§) is similar to the t — J model. 
Noting that t'/so is a small expansion parameter, we can 
regard Eqs. (|J) and (||) as a perturbation to Eqs. (||) 
and (Q) ; the latter parts are essentially a spinless fermion 
Hamiltonian with nearest-neighbor hopping. The degen- 
eracy with respect to spin can be lifted only after Eqs. (||) 
and (^) are introduced. 

As far as t' /eq is small, we can apply the perturba- 
tion treatment which was used for the large-?7 limit of 
the Hubbard model. |2(J Just like in the latter case the 
ground state wave function |^ s ) is given as a direct prod- 
uct of spinless fermion wave function |$sf) and spin wave 
function \x)'- 



|* 9 > = |$sf)®|x>- 



(37) 



To discuss the ground state of the effective Hamiltonian 
we take for |$sf) the spinless fermion ground state of 



Ti-SF = s}Xa]a j+1 + H.c), 



(38) 



where aj and aj are spinless fermion annihilation and cre- 
ation operators respectively. |x) has to be determined to 
minimize the spin-dependent part of the effective Hamil- 
tonian. Taking the average of 7i e ff over |$sf) and picking 
up the spin-dependent part, we obtain 



(^eff)sF — "/cff^S,; • Sj+i, 

i 



where 



1 - 


/sin(7w)\ 2 " 




V TTV ) . 



g{v) = v 1 



f{v) = — smlitv) ttu cos(ifv) — sin(7w) 



h{v) 



v 2 [\ 



sin(27w) 
2-kv 



(39) 
(40) 

(41) 
(42) 
(43) 



The summation in Eq. (39) is taken over the squeezed 
spins. |2(J v is twice the density n of electrons. For 
< v < 1, g[v) and h(y) are positive, while j[v) is 
negative. Therefore, for positive s and t, a competition 
occurs among three terms: the s(t' /so) 2 term favors fer- 
romagnetism, whereas the t(t' /sq) 2 and s 2 /U terms tend 
to destabilize it. Fig. |ll| shows the sign of J c ff. For 
< v < 1, the ground state is ferromagnetic as long as 
Jeff is negative. 

Now let us turn to the excitation spectrum of the 
one spin flipped state in the same parameter region of 
Jeff "C t. The low energy excitation of the spin part are 
the spin waves with momentum q = 2ttJ/N & and energy 
e(q) = J c ff(l — cosg). The momentum of the spin part 
changes the boundary conditions of the spinless fermions, 
as noted by Woynarovich, p^] so that the momenta of 
spinless fermions are (L/2)kj = 2nlj +q with Ij integers, 
and this phase shift appears in the energy and momen- 
tum of the spinless fermion part: 



E = J e ff(l — cosq) — 2s cos 



2ttL 



(L/2) 



P = ^k j =2n-J2l: 



2N K 
L 



-q. 



(44) 



In the thermodynamic limit, if the ground state of q = 
(i.e. of the highest spin state) is nondegenerate, then the 
energy E{P) — E-pu of the lowest lying excitations of the 
system is given by 



oj(k) = J c ff(l — cosg) H sq 2 



(45) 



where k = vq+PpM- Unlike the insulating case at n = i, 
this spin wave excitation is not separated any more by a 
finite gap from the Stoner continuum. We do not know 
the relevance of this fact to the stability of itinerant fer- 
romagnetism in this case. 

Depending on the number of particles, see Table. 0, 
it can easily happen that the fully polarized state is de- 
generate. In this case of orbital degeneracy the energy of 
the spin wave with momentum k — vq ± |Pfm| becomes 

2 , x0 sin7ri/ 
u(k) = J eff (l-cosg) + -s(g=F7r) 2 (46) 

L 7T 

where we can clearly identify the competition between 
the kinetic energy gained by twisting the boundary con- 
ditions and the energy loss due to the spin wave. Ac- 
tually, in the latter case, unless J ff > J c oc s^sin7ri/ 
, the minimum energy will occur for some finite q, and 
it prevents the highest spin state from being the ground 
state. 

We have performed an exact diagonalization study of 
the model at n = 1/4 (i.e. v = 1/2) to determine the re- 
gion where the fully polarized is stable against spin flip. 
The results are presented on Fig. O. The behavior de- 
scribed above can be clearly seen for larger values of U 
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and |eo|, whereas near the flat band for £o > £fb the 
effect of orbital degeneracy is not so important . 



rahmanyam and X. Zotos. We also thank IDRIS (Orsay) 
for allocation of CPU time on the C94 and C98 CRAY 
supercomputers. 



ACKNOWLEDGMENTS 

We are especially thankful to H. Tasaki for stimulat- 
ing discussions. The authors enjoyed interesting discus- 
sions with P. Fazekas, K. Hallberg, J. Solyom, V. Sub- 



APPENDIX: HIGHER ORDER TERMS IN THE 
SERIES EXPANSION 

Expansion of K\ continues as: 



id 4) = - 



16s 4 At' 4 



8t' 



U 3 ep e§(2e o + C0' 



id 5) - - 



12s 3 i' 2 40st' 4 64s 3 t' 2 32s 3 t' 2 16st' 4 



36s 2 i' 4 120s 6 



s U 3 
96s¥ 4 



ep 2 

4j./2 



+ 



kP = 



16s 2 t 



2 + /4 



■ 2 u 3 

32sH' 4 



64s 4 i' 
~ep 3 



ep 

4m/2 



64s 4 t 



2 + /4 



elU 2 



96s 2 t 

'du 2 ' 



56sH' 2 

4u 



16s 2 t' 4 

ep 



A8sH 



2 + /4 



96t 



e 2 (2e + U f 
84? 720F 



eg (2e + Uf 
360F 136F 



soU 5 



ep 4 



ep 3 



28? 



32? 



ei (2e 
236F 

80? 



- U) e 4 (2e + U) ' 
224? 



4u 



32? 



4(e + U) e 3 (2e + U) 3 e 4 (2e + U) 2 e% (2e + U) 
where t = t' = s. Contribution to K 2 starts in 4th order: 



16t' 6 



+ 



8sH 



4 + /2 



et (so + U) 



4s 4 
U 3 

16s 3 *' 2 
s U 3 
8sH' 4 



8sH 



3 + /2 



e 2 U 2 ' 



48s 6 16sV 2 24sV 4 16s 4 t' 2 24s 2 *' 4 20s 4 t' 2 4s 2 t' 4 



K. 



(7) 



288F 



+ ■ 



u 5 
iu? 



P 3 



P 3 



P 2 



P 2 



16? 32? 96? 



32? 



e U 5 e U 4 



ep 3 



ep 2 



ep 



+ 



ep 
16? 



ep 



e 5 (s + U) el (2e + U) ' 



4sH 



4-1/2 



AsH'" 



e 4 (e + U) e 4 (2e + U) ' 



and to K3 in 6th order: 
3 ~W 

(7) _ 48F 24F 8F 12F 
3 ~ ~e~olP ~ Tp 4 ~ ep 3 ~ ep 2 ' 

The coefficients Kj for U, — £0 3> t, where t = s = t' 
and we introduced the notation £0 = — £0, are: 



(6) 



20F 24i 6 



16t b 



i6r 



?P el(2e Q + U) 2 s 4 (2e- a + Uy 



(7) _ 72F ^ 32/ 



?7 



160F 



352F 



f6 



(2£ + Uf 4 (2£ + U) ' 



and the first contributions to Kq come in 6th order: 
R{6) 4f 6 16t 6 

A 2 - + 



£ 4 (7 £ 4 (2£ + U) 1 



fc 



8F 

£"5[/ 



48F 



48F 



£ 4 (2£ + C/) 2 el (2£o + C/) ' 
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TABLE I. £fb, constant C and critical values of Coulomb repulsion U for some selected values 
oit'/s. 



p > p < 



t'/s 


efb/s 


C 


Uo/p 


Uc/P 


U* pp /p 


Ud/\ P \ 


Uc/\p\ 


U™/\p\ 


4 


14 


0.00154 


0.830 


0.831 


0.223 


2.402 


2.414 


0.271 


2 


2 


0.0136 


1.047 


1.077 


0.707 


1.851 


1.945 


1 


V2 





0.0302 


1.235 


1.355 


1.155 


1.527 


1.715 


1.406 


1 


-1 


0.0557 


1.476 


1.860 


1.789 


1.277 


1.555 


1.505 


1/2 


-1.75 


0.149 


2.169 


3.88 


3.88 


1.025 


1.294 


1.294 



TABLE II. Momentum Pfm (with degeneracy) of the fully 
polarized state for periodic (PBC) and antiperiodic (APBC) 
boundary conditions, t — 0. 



£() > £FB £0 < £FB 



L/2 


N e 


PBC 


APBC 


PBC 


APBC 


even 


even 


±vn 










even 


odd 


TV 


±(1-z/)tt 







odd 


even 





±!/7T 







odd 


odd 


±(1-v)tt 


7T 





±VTX 
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FIG. 1. Schematic representation on the ID model. Solid 
circles stand for 'c' sites, and the open ones for 'd' sites. 
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FIG. 2. Some possible realizations of the model in two 
dimensions: 'c' sites form square (a), honeycomb (b) and 
kagome (c) lattice. 



0.01 



-0.01 



-0.02 



-0.03 




FIG. 5. The different Pade approximants for K\ and 
K2 are shown for U/eo = 8. The series is convergent for 
s/so ~ 0.15. The dashed line represents K[ 2 ^ + K[ 3 \ Note 
the excellent agreement with the values extracted from ex- 
act diagonalization data of a 16 site cluster (diamonds and 
triangles). The finite size effects are negligible. 
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FIG. 3. Phase diagram for t = 0, s = t' = 1 and quarter 
filling. A flat band is realized for £0 — — 1 an d for U = +00 
the model is ferromagnetic for any eo- For U < U c the system 
is paramagnetic, with strong antiferromagnetic correlations. 
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FIG. 6. Stability of the ferromagnetic state for s = t' 
and t = 0, quarter filling. Solid lines for large values of U c 
come from Pade approximants of the perturbational expan- 
sion (main figure* and insert). Solid straight lines are the 
approximation near the flat band. Opens symbols are the 
data for local stability for system sizes 12, 16 and 18 with pe- 
riodic boundary condition. The finite size effect is negligible 
apart from the effect of degeneracy near e = —2s for 12 and 
16 site systems. Solid points are global stability data of the 
16 site cluster. 



FIG. 4. The third order process leading to ferromagnetic 
effective exchange J e g — Ast' 2 /eq. 
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FIG. 7. The excitation spectra in the thermodynamic 
limit for some selected cases, (a) For small values of U the 
fully polarized state is not stable, but already we can see a 
well defined magnon, which is very close to the particle-hole 
excitation continuum for q = it. (b) For larger values of U we 
can see a normal ferromagnetic magnon dispersion relation, 
with the Stoner continuum pushed up to higher energies. 
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FIG. 8. The comparison of the excitation spectrum ob- 
tained from matrix M p _ p / (q) (solid diamonds) to exact di- 
agolization of a 16 site cluster (open squares). Although 
U/s — 0.3 is already comparable to the band gap, the errors 
are very small. The discrepancy is due interband scattering 
and it scales as U 2 . Here k = ir — q, where it is the momentum 
of|FP). 
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FIG. 9. U c /\p\ calculated from exact diagonalization 
(open symbols) and from diagonalization of matrix M PiP /(7r) 
(solid symbols) for system sizes L — 12, 14, 16 and 18 with 
periodic boundary conditions, t' = s and t — 0. The solid 
circle stands for the value of U c in the thermodynamic limit 
L — > oo. 
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FIG. 10. (a) The dispersion of the lowest energy spin 
nipped state for t' = s = 1, £o = — 1 and t = from exact di- 
agonalization of 12,16, and 20 site clusters with periodic (P) 
and L — 14 and 18 with antiperiodic (A) boundary condi- 
tions. Although the S max state at k — it is higher in energy, 
the dispersion has clearly ferromagnetic k dependence, (b) 
The finite size scaling of the energy difference, Egs is the 
lowest energy in the S z = Smax — 1 subspace. 
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FIG. 11. For U = +oo and t'/e < 1 the sign of J eS 
and the nature of the ground state is determined by the ratio 
t/s. For negative t it will always give a ferromagnetic effective 
coupling. 
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FIG. 12. Local stability for n = 1/4 from exact diagonal- 
ization of 12, 16, 20 and 24 site clusters with periodic (P) or 
antiperiodic (A) boundary conditions (t' = s, t = 0). The 
fully polarized state is degenerate for open symbols. Ferro- 
magnetism is stable against single spin flip for U > U c . 
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